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We investigate the possibility to determine the fraction 
of compact objects in the Universe by studying gravitational 
lensing effects on Type la supernova observations. Using sim- 
ulated data sets from one year of operation of the proposed 
dedicated supernova detection satellite SNAP, we find that 
it should be possible to determine the fraction of compact 
objects to an accuracy of < 5%. 



I. INTRODUCTION 

Recent measurements of anisotropics in the cosmic mi- 
crowave background radiation (CMBR) show that the 
Universe is very close to fiat, i.e., f2 to t ~ 1 (jjj^) ■ Since ob- 
servations of Type la supernovae (SNe) indicate that the 
expansion rate of the Universe is accelerating, the major 
part of this total energy should have negative pressure, 
e.g., in the form of the cosmological constant correspond- 
ing to Qa ~ 0.7 in a flat universe ||,||], in agreement 
with constraints on the matter density J7m from cluster 
abundances |HJ|(| and large-scale structure 0. Thus, a 
a concordance model with J7m ~ 0.3 and J7a ~ 0.7 has 
emerged. 

The constitution of the total matter density is a mat- 
ter of intense theoretical and experimental research. The 
energy density in baryonic matter as derived from Big 
Bang nucleosynthesis (BBN) is given by Q,\,h 2 = 0.019 ± 
0.0024 H whereas CMBR measurements yield a some- 
what higher value; 0.022 < tt h h 2 < 0.040 [§. Either 
way, the matter density in baryons is almost one order 
of magnitude smaller than the non-baryonic dark mat- 
ter (DM) component. However, the BBN range for the 
baryon density still means that most of the baryonic mat- 
ter is also dark. 

There are various possibilities for the dark baryons 
to hide, examples are warm gas in groups and clusters 
(which is difficult to detect at present jT(i|j), or in the 
form of massive compact halo objects (MACHOs), where 
indeed there have been detections p] , ^2[ . The long lines 
of sight to distant supernovae means that they are well 
suited to probe the matter content along the paths of the 
light rays, where gravitational lensing may occur where 
there are matter accumulations. Compact objects give 
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more distinct lensing effects, enabling a distinction be- 
tween diffuse matter and compact bodies along the light 
path. In particular, it may be possible to investigate 
whether the halo fraction deduced for the Milky Way 
from microlensing along the line of sight to the Large 
Magellanic Cloud, of the order of 20 % [[l), is a uni- 
versal number or if the average cosmological fraction is 
larger or smaller. 

Regardless of its constitution, we can classify dark 
matter according to its clustering properties. In this pa- 
per we will use the terminology of smooth DM for DM 
candidates which tend to be smooth on subgalaxy scales, 
e.g., weakly interacting massive particles (WIMPs) such 
as neutralinos. The term compact DM will be reserved for 
MACHOs such as brown or white dwarfs and primordial 
black holes (PBHs) §J|. 

An advantage with gravitational lensing is that its ef- 
fects can determine the distribution of dark matter inde- 
pendent of its constitution or its dynamical state. The 
topic of this paper concerns the use of the gravitational 
magnification of standard candles, such as Type la SNe, 
to determine the fraction of compact objects in a cosmo- 
logical context. 

In an early study, Rauch jl4|] concluded that with a 
sample of 1000 Type la SNe at redshift one should 

be able to discriminate between the extreme cases of all 
DM as smooth or in compact objects. More recent work 
(Metcalf and Silk Goliath and Mortsell has 
shown that a sample of 50-100 SNe should be sufficient to 
make the same discrimination. Seljak and Holz [fl7| have 
found that it should be possible to actually determine 
the fraction of compact objects to 20 % accuracy with 
100-400 Type la SNe at z = 1. In this work we extend 
these studies by exploring the possibility to determine 
the fraction of compact objects using a future sample of 
Type la SNe distributed over a broad redshift range, with 
the intrinsic spread of absolute luminosity of the SNe and 
expected measurement error for a proposed space-borne 
mission, SNAP [pi taken into account. 

In Sec. |n] we clarify the distinction between compact 
and non-compact dark matter objects. In Sec. [Hi] we 
discuss gravitational lensing of SNe and in Sec. |IV| , we 
present the Monte-Carlo simulation package used to pre- 
dict this effect. Sec. |v| is concerned with the proposed 
satellite telescope which will be able to produce the data 
sets used in this study and in Sec. VI we present our 
results. The paper is concluded with a summary in 
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II. COMPACT OBJECTS 

For an object to be compact in a lensing context, we 
demand that it is contained within its own Einstein ra- 
dius, ve- For a DM clump at z — 0.5 and a source at 
z= 1, r B ~ lO^M/M©) 1 / 2 pc in a fl M = 0.3,0a = 0.7 
and h — 0.65 cosmology. Also, the Einstein radius pro- 
jected on the source plane should be larger than the size 
of the source. A Type la SN at maximum luminosity 
has a size of ~ 10 15 cm, implying a lower mass limit of 
~ 10 _4 M Q for the case described above. 

The effects of lensing by compact objects arc different 
from those of lensing by halos consisting of smoothly dis- 
tributed dark matter, such as in the singular isothermal 
sphere or Navarro-Frenk- White density profile (NFW; 
p9[), the main difference being the tail of large mag- 
nifications caused by small impact-parameter lines of 
sight near the compact objects. However, N-body sim- 
ulations predict, besides the overall cuspy profiles of 
ordinary Galaxy-sized dark matter halos, also a large 
number of small subhalos on all length scales which 
can been resolved p9| , p0| . The number density of the 
smaller objects, of mass M, follows approximately the 
law dN/dM oc M~ 2 as predicted by Press-Schechter the- 
ory. Thus one may expect a multitude of subhalos in 
each galaxy or cluster halo. In addition, N-body simu- 
lations show the less massive halos to be denser (mainly 
due to them being formed early when the background 
density was higher). Thus, it is appropriate to address 
the question whether this type of small-scale structure, 
and in particular the dense central regions of them, may 
give rise to lensing effects similar to the ones caused by 
truly compact objects. 

To put a bound of these possible effects, we use the re- 
sults of the most accurate numerical simulations to date 
pot . The density profiles within dark matter clumps ob- 
tained in the simulations can be fit by the Moore profile 
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Here p' M and a are not independent parameters, but 
related by the concentration parameter cm = R200/&, 
which depends on mass roughly as 



cm 



10.6Mr° 084 . 



(2) 



(Here Mi 2 is the virial mass in units of 10 M©; -R200 
is the virial radius where the average overdensity is 200 
times the background density. Similar relations appear 
in the NFW simulations.) 

Using these relations we can derive the mass Mr (also 
in units of 10 12 M©) within distance R pc parsecs from 
the center of the halo, 



M R ~ 10- 7 M° 2 4 -Rp C 5 



(3) 

Comparing this with the Einstein radius for the same 
mass for a lensing event of typical distance Dq pc Giga- 
parsecs, 



R E ~ 10 4 M° 5 ^ C pc, (4) 
we find that Mr is within its own Einstein radius if 



Mr < 10 



m 12 u Gpc- 



(5) 



The requirement that the lensing mass is greater than 
10~ 4 M© means that dark matter clumps of the Moore 
type with mass greater than around 10 4 M© (A/12 > 
10~ 8 ) will in principle contribute to compact lensing. 
However, it is only the very central, dense core that can 
contribute to the lensing, and we see from Eq. (g) that 
only a small fraction / of the mass in the clump is within 
its own Einstein radius, 
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lQ- 4 M™D 3 Gpc . 



(6) 



Thus only clumps of Galaxy size contribute appreciably 
to the lensing, and this part is included in our standard 
calculations. 

Since this analysis was made for Moore- type halos, 
which are more concentrated than NFW halos, we con- 
clude that compact objects detected through lensing of 
supernovae cannot, according to current thinking about 
structure formation, be caused by clumps of particle dark 
matter formed through hierarchical clustering. 



III. GRAVITATIONAL LENSING OF 
SUPERNOVAE 

The effect of gravitational lensing on Type la SN mea- 
surements is to cause a dispersion in the Hubble diagram. 
In Fig. ^ we compare the dispersion due to gravitational 
lensing with the intrinsic dispersion and the typical mea- 
surement error for Type la SNe. In the upper left panel, 
we show an ideal Hubble diagram with no dispersion and 
in the upper right panel we have added the dispersion due 
to lensing (lens) in a universe with 20 % compact objects 
and 80 % smooth dark matter halos parametrized by the 
Navarro-Frenk- White formula ]l9fl . Comparing with the 
panel in the lower left where the intrinsic dispersion (intr) 
and measurement error (err) have been included, we see 
that the effects become comparable at a redshift of unity. 
In the lower right we see the most realistic simulation 
with intrinsic dispersion, measurement error and lensing 
dispersion. 
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FIG. 1. A comparison of the dispersion due to lensing (up- 
per right) and the measurement error and intrinsic dispersion 
of Type la SNe (lower left). 

Of course, the additional dispersion caused by gravi- 
tational lensing will be a source of systematic error in 
the cosmological parameter determination with Type la 
supernovae. However, a possible virtue of lensing is that 
the distribution of luminosities might be used to obtain 
some information on the matter distribution in the Uni- 
verse, e.g., to determine the fraction of compact DM in 
our Universe. 



The exact parameterization of the smooth DM halo pro- 
file is not important for the results obtained in this paper 
p2fl . The results are also independent of the individ- 
ual masses of the compact objects as well as their large 
scale clustering properties [^l|,^2). As we have seen, the 
eventual small-scale structure in the "smooth" compo- 
nent does not act as a compact component. 

V. SUPERNOVA/ ACCELERATION PROBE 

To make realistic predictions of the statistics and qual- 
ity of the supernova sample, we use the projected dis- 
covery potential of the Supernova/ Acceleration Probe 
(SNAP) project. This is a proposed satellite telescope 
capable of discovering over 2000 T ype la SNe per year in 
the redshift range 0.1 < z < 1.7 jig] . The most antici- 
pated use of the data is to gain further accuracy in the 
determination of, e.g., fijvi and the curvature of the Uni- 
verse fik, but also to give insight into the nature of the 
negative pressure energy component by constraining the 
equation of state or the redshift dependence of 

the effective energy density [p5| . Here we show how the 
data also can be used to give information on the fraction 
of compact objects in the dark matter component. 

More specifically, it is projected that in one year, 
SNAP will be able to discover, follow the light curves and 
obtain spectra for of the order of 2000 SNe. While the 
exact redshift distribution of the Type la supernovae to 
be followed by SNAP might be changed upon studies of 
the optimal search strategies for the primary goals of the 
project, we have used the Monte-Carlo generated sample 
spectra of 2366 SNe distributed according to the SNAP 
proposal ]l8| , see Fig. to make specific predictions. 



IV. SUPERNOVA OBSERVATION CALCULATOR 

To perform realistic calculations we have developed 
a numerical simulation package, SNOC (the SuperNova 
Observation Calculator). It can be used to estimate var- 
ious systematic effects such as dust extinction and grav- 
itational lensing on current SN measurements as well as 
the accuracy to which various parameters can be mea- 
sured with future SN searches. In this paper, we use 
SNOC to obtain simulated samples of the intrinsic dis- 
persion and gravitational lensing effects of Type la SNe 
over a broad redshift range. 

The intrinsic dispersion and measurement error is rep- 
resented by a Gaussian distribution with a m = 0.16 mag. 
Gravitational lensing effects are calculated by tracing the 
light between the source and the observer by sending it 
through a series of spherical cells in which the dark mat- 
ter distribution can be specified. For more details on 
the method, originally proposed by Holz and Wald, see 

HE§. 

We will model compact DM as point-masses and 
smooth DM as the Navarro-Frenk- White density profile. 
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FIG. 2. Number of Type la SNe expected for various red- 
shift bins in a one-year exposure with the proposed SNAP 
satellite. The data are taken from Table 7.2 in the SNAP 
science proposal 



VI. RESULTS 

Using SNOC, we have created large data sets of syn- 
thetic SNe observations with a variable fraction of com- 



3 



pact objects ranging from to 40 % using the follow- 
ing cosmological background parameter values: Qm = 
0.3, f^A = 0.7 and h = 0.65. (For a discussion of how the 
halo distributions were generated, see |22].) These are 
used as reference samples. In Fig. |^, we have plotted the 
dispersion in the reference samples for % (full line), 20 
% (dashed line) and 40 % (dotted line) compact objects. 
In the upper panel, the lensing dispersion is plotted, i.e., 
the dashed line basically corresponds to the scatter in the 
upper right panel in Fig. [l| The zero value corresponds 
to the value one would obtain in a homogeneous universe 
(upper left panel in Fig. Q). Note that negative values 
correspond to magnified events, positive values to demag- 
nified events. It is clear, that as the fraction of compact 
DM grows, lensing effects becomes larger in the sense 
that we have a larger fraction of highly magnified events. 
In the lower panel we have added a Gaussian intrinsic dis- 
persion and measurement error, cr m = 0.16 mag, making 
the distributions look more similar (cf., lower right panel 
in Fig. 0). Although this smearing obviously decreases 
the significance of the compact signal, it can be seen that 
the high-magnification tails still are clearly visible. 



mental sample with our high-statistics reference samples 
using the Kolmogorov-Smirnov (KS) test, we obtain a 
confidence level for the hypothesis that the experimental 
sample is drawn from the same distribution as the ref- 
erence sample. For each fraction of compact objects in 
the experiments (6, 11 and 21 %), we repeat this proce- 
dure for 1000 experimental realizations and pick out the 
reference sample which gives the highest confidence level 
for each experiment. Plotting the number of best-fit ref- 
erence samples as a function of the fraction of compact 
objects in the reference sample, we can fit a Gaussian 
and thus estimate the true value and the dispersion. In 
each case, we get a mean value within 1 % of the true 
value and a one-sigma error less than 5 %, see Fig. 0. 
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FIG. 3. Magnitude dispersion of reference samples for % 
(full line), 20 % (dashed line) and 40 % (dotted line) compact 
objects. The bottom panel includes a Gaussian smearing, 
<7 m = 0.16 mag, due to intrinsic brightness differences between 
supernovae and from measurement error. The distributions 
show the projected scatter around the ideal Hubble diagram 
for Type la supernovae with a relative redshift distribution as 
in Fig. |. 

We have also created a large number of simulated one- 
year SNAP data sets (according to Fig. || above) with 
6, 11 and 21 % compact objects. These are our exper- 
imental samples. By comparing each generated experi- 
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FIG. 4. The number of best-fit reference samples as a func- 
tion of the fraction of compact objects in the reference sample. 

As one could expect, and is shown in Fig. [j], lensing 
effects get larger at higher redshifts. At low redshifts 
the dispersion is completely dominated by the intrinsic 
dispersion. Therefore, we have only used the data from 
SNe at z > 0.8, a total of 1387 SNe, to obtain the result 
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in Fig. |Q 

In this analysis, we have assumed that Om and fl\ 
will be known to an accuracy where the error in luminos- 
ity is negligible in comparison to the intrinsic and lcnsing 
dispersion of Type la SNe. This assumption is not unrea- 
sonable with future CMBR observations combined with 
other cosmological tests and the SNAP data itself, nor is 
it crucial in the sense that we are dealing with the disper- 
sion of luminosities around the true mean value, not the 
mean value itself. In order to test the sensitivity of our 
results to changes in the cosmological parameter values, 
we have performed a number of Monte-Carlo simulations 
using different sets of parameters in our experimental and 
reference samples and found the error in the total energy 
density in compact objects to be negligible^. Also, we 
have used a value of the intrinsic dispersion and mea- 
surement error of a m = 0.16 mag, a value which may 
be appreciably smaller in the future when the large data 
sample expected may allow, e.g., a more refined descrip- 
tion and classification of supernovae. Using simulated 
data sets with a m = 0.1 mag, the one-sigma error in the 
determination of the fraction of compact objects, using 
the same cosmology as above, becomes less than 3 %, as 
depicted in Fig. [|. 



VII. SUMMARY 

The proposed SNAP satellite will be able to detect and 
obtain spectra for more than 2000 Typa la SNe per year. 
In this paper we have used simulated data sets obtained 
with the SNOC simulation package to show how one- 
year SNAP data can be used to determine the fraction 
of compact DM in our Universe to < 5% accuracy, as- 
suming the intrinsic dispersion and measurement error is 
(T m = 0.16 mag. If the intrinsic dispersion and measure- 
ment error can be further reduced, e.g., from a improved 
understanding of Type la SN detonation mechanisms, the 
accuracy can be improved even further. 
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FIG. 5. The number of best-fit reference samples as a func- 
tion of the fraction of compact objects in the reference sample 
for an intrinsic dispersion and measurement error of <r m = 0.1 
mag. 



1 Of course, the data from SNe at z < 0.8 can still be used to 
constrain the values of £Im and Qa- In fact, with one year of 
SNAP data, it is possible to determine Qm with a statistical 
uncertainty of Af2 M ~ 0.02 @. 

2 Note that the effect from lensing is proportional to the total 
energy density in compact objects, not the fraction of compact 
objects, see, e.g., |^6j. A higher Qm would therefore increase 
the sensitivity for smaller fractions. 
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